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(57) Abstract 

The invention relates to a communication system (300) having a receiver (304) that is capable of performing targeted interference 
suppression. An interference clasisifier (314) within the receiver (304) analyzes a signal received from a channel (306) and identifies and 
classifies interference components within the signal. An interference suppressor (316) then suppresses the interference components in the 
signal based on interference type. The interference suppressor (316) includes a plurality of interference suppression modules that are each 
optimal for suppressing certain interference types. In one embodiment, one or more of the interference suppression modules perform 
interference suppression using modal moment estimates to identify different modes within the signal. 
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INTERFERENCE SUPPRESSION USING MODAL MOMENT ESTIMATES 



Field of the Invention 



10 



The invention relates in general to communications systems and, 
more particularly, to communications systems having Interference mitigation 
capabilities. 

Backgroun d of the Invantinn 



Digital cellular communications systems, paging systems, land mobile 
radio, and mobile battlefield communication systems are called upon to 
operate effectively in increasingly adverse spectral environments. Obstacles 
to digital mobile communications include non-stationary co-site Interference, 

15 ambient communication signals, and hostile Jamming Interference. Typical 
interference levels can be very high, especially in a tactical environment, 
and can often exceed the desired signal by 60 decibels or more. Such 
levels ovenvhelm the operational capabilities of current tactical and 
commercial radio systems. Furthermore, PCS, paging, and cellular 

20 communication networks that Incorporate spread-spectrum technologies, 
such as CDMA, increasingly face capacity limitation caused in part by 
growing spectral clutter. The presence of ambient high power partialband 
and narrowband signals often precludes frequency reuse by spread- 
spectrum communication systems. 

25 Therefore, there is a need for a method and apparatus for reducing 

thfe impact of undesirable spectral components in a communications system. 

Brief De-sn ription of the Drawings 

f^'G- 1 's a high level block diagram illustrating a communication 
system in accordance with one embodiment of the present invention; 
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FIG. 2 is a block diagram illustrating a communication system In 
accordance with another aspect of the present Invention; 

FIG. 3 is a block diagram illustrating elements within the disturbance 
classification unit 74 of FIG. 2; 
5 FIG. 4 is a joint time frequency (JTF) feature plane diagram illustrating 

Inbound communication signal interferers versus time; 

FIG. 5 is a cross section of a feature plane transfomiation diagram for 
a five FSK interference signal; 

FIG. 6 is a flowchart Illustrating a transfomiation masking method in 
1 0 accordance with one embodiment of the present invention; 

FIGS 7-1 1 are graphs Illustrating a method for performing 
preclasslfication filtering in accordance with one embodiment of the present 
Invention; 

FIGS. 12 and 13 are cross-sections of feature plane signal 
15 transformations illustrating the beneficial effects of using phase domain 
filtering; 

FIG. 14 is a spectrum diagram illustrating a method for detemiining a 
baud rate for a binary phase shift keying (BPSK) signal in accordance with 
one embodiment of the present Invention; 
20 FIG. 15 illustrates a linear array of alternating Kronecker deltas that Is 

used to determine a baud rate for an frequency shift keying (FSK) signal in 
accordance with one embodiment of the present invention; 

FIG. 16 is a spectrum diagram illustrating a method for determining a 
baud rate for an FSK signal in accordance with one embodiment of the 
25 present invention; 

FIG. 17 is a plot illustrating modal moment estimates for a data signal 
corrupted by a partial band noise jammer in accordance with one 
embodiment of the present invention; 

FIG. 18 is a flowchart illustrating a process for deriving modal moment 
30 estimates In accordance with one embodiment of the present invention; 
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FIG. 19 is a spectrum diagram illustrating a r)umber of bands WsN for 
use In frequency hopping In accordance with one embodiment of the 
present Invention; 

FIG. 20 is a flowchart illustrating a method for performing frequency 
5 hop masking In accordance with one embodiment of the present inventlon; 

FIG. 21 Is a flowchart Illustrating the operation of an inverse whitening 
function In accordance with one embodiment of the present invention; 

FIG. 22 Is a spectrum diagram illustrating the spectrum of a spread 
spectrum data signal corrupted by a PSK partial band jammer 122; 
"•0 FIG. 23 is a flowchart Illustrating the operation of an adaptive inverse 

weight interference suppression method In accordance with one 
embodiment of the present Invention; 

FIG. 24 is a flowchart Illustrating the operation of an adaptive excision 
Interference suppression method in accordance with one embodiment of the 
15 present invention; 

FIG. 25 is a flowchart Illustrating a cascaded Interference suppression 
procedure in accordance with one embodiment of the present invention; and 

FIGS. 26-28 are spectrum diagrams Illustrating a cascaded 
suppression procedure in accordance with one embodiment of the present 
20 Invention. 

Detailed Description of the Drawings 

The present invention relates to a communications system that Is 
25 capable of reducing the impact of undesirable spectral components (e.g., 
noise, Interference, and jamming) In a receive signal. The system includes 
an Interference classifier for classifying each of the undesired spectral 
components (if any) by interference-type and an interference suppressor for 
suppressing the undesired interference components In the receive signal 
30 using one or more interference suppression methods that may be "targeted" 
to the specific type of Interference. By targeting Interference suppression by 
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interference-type, a superior level of interference suppression performance 
Is achieved as compared to prior art systems that use a common 
suppression method for all types of interference. In one embodiment of the 
invention, the interference classification and suppression functions are 
5 performed automatically. 

The Invention is capable of operating in real time on a short sample 
set, with power consumption that is within the capabilities of modern digital 
signal processors (DSP). In addition, the invention is capable of optimizing 
bit error rate (BER) performance given a variety of interferers without 
1 0 sacrifice of bandwidth or system capacity and, in some cases, can improve 
the bit error rate (BER) of a communications system by multiple orders of 
magnitude. The invention is of relatively low complexity and can be 
implemented using computationally efficient recognition/cancellation 
routines. In one embodiment, the invention is used to provide for frequency 
1 5 reuse in a situation where a spread spectrum communications system Is 
overiaying one or more narrowband and/or partialband systems. 

As described In more detail below, the principles of the present 
invention can be, and preferably are, implemented in a spread spectrum 
communications system. A spread spectrum communications system is a 
20 system where the bandwidth of the transmitted radio frequency signal is 
wider than that required by the data rate and modulation type of the 
underlying Infomiation signal. That is, a second layer of modulation Is used 
that "spreads" the spectrum of the modulated Information signal to provide a 
number of important advantages. For example, spread spectmm systems 
25 are generally better at rejecting interference than are other types of 

communications systems. In addition, spread spectrum systems can support 
the use of code division multiple access (CDMA) techniques to provide a 
plurality of separate communication channels In a given bandwidth. Another 
advantage of spread spectrum techniques is that they are conducive to 
30 signal hiding in secure communications applications. Spread spectrum 

systems are also capable of high resolution ranging. Various methods, such 
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as direct sequence spread spectrum (DSSS) and frequency hopping spread 
spectrum (FHSS), can be used to spread the bandwidth of the transmit 
signal. 

FIG. 1 is a high level block diagram illustrating a system 300 in 
5 accordance with one embodiment of the present invention. As illustrated, 
the system 300 Includes a transmitter 302 which communicates with a 
receiver 304 through a channel 306. The transmitter 302 receives data from 
a data source 308 and processes the data to create a transmit signal that is 
delivered to the channel 306. The receiver 304 receives a signal from the 

10 channel (which is a modified fomn of the transmit signal) and processes the 
signal to recover the original data. The recovered data Is then delivered to 
the data sink 310. 

The transmitter 302 includes, among other things, a 
modulation/coding unit 312 that applies modulation and/or coding to the 

15 data from the data source 308. For example, the modulation /coding unit 
312 can apply source coding, channel coding, interleaving, and/or 
upconverslon to the data signal. In a spread spectrum system, the 
modulation/coding unit can also apply signal spreading to the signal using 
methods that are well known in the art. 

20 The receiver 304 includes an interference classifier 314, an 

interference suppressor 316, and a demodulation/ decoding unit 318. The 
receiver 304 receives the signal from the channel 306 in a signal receptor 
(not shown), such as an antenna. The interference classifier 314 analyzes 
the signal received from the channel 306 and Identifies and classifies 

25 interference components within the signal. The interference components 
can be from any of a number of different sources, such as nearby 
communications systems and/or hostile entities attempting to jam 
transmissions from the transmitter 312. The interference classifier 314 
outputs a signal indicative of the interference classification of each of the 

30 identified interference components. The interference suppressor 316 
receives both the classification signal from the interference classifier 314 
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and the signal received from tlie channel 306. The interference suppressor 
316 then selects one or more suppression methods that work best for the 
types of interference Identified and uses the selected methods to suppress 
interference in the receive signal. After the interference components have 
5 been suppressed, the receive signal is transferred to the 

demodulation/decoding unit 318 to remove all modulation and/or coding 
from the signal. The signal is then delivered to the data sink 310. 

FIG. 2 is a more detailed block diagram Illustrating a spread spectrum 
communications system 10 in accordance with one embodiment of the 
0 present invention. As illustrated, the system 10 includes an input device 12, 
a transmitter 14, a channel 16, a receiver 18, and an output device 20. The 
transmitter 14 receives an Input signal from the input device 12 and 
processes the signal into a form for transmission into the channel 16. The 
channel 16 carries the signal from the transmitter 14 to the receiver 18. 
5 Receiver 18 processes the signal received from the channel to recover the 
original information that was generated by the input device 12. This 
recovered infomiation is then transferred to the output device 20. 

The input device 12 can Include virtually any type of information 
source. That is, the input device 12 can provide, for example, audio 
) Infomiation (such as speech), computer data, or video infonnation. In the 
illustrated embodiment, input device 12 comprises a microphone for 
converting speech Into an analog electrical signal indicative thereof. 

■ The transmitter 14 includes: analog to digital (A/D) converter 22, 
encoder 24, multiplier 26, digitally controlled oscillator (DCO) 28, N level 
' chip sequence generator 30, chip shaping filter 32, frequency hopping (FH) 
adaptation module 34, and processing gain adaptation module 36. A/D 
converter 22 is operative for digitizing the analog signal from the Input 
device. That Is, the A/D converter 22 samples the analog signal from the 
input device 12 and creates digital values (i.e., samples) that are each 
Indicative of the magnitude of the analog signal at a particular instant in time 
(I.e., the sampling time). It should be appreciated that an A/D converter 22 is 
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only necessary if the Input device 12 Is providing analog Input. The encoder 
24 receives the digital signal from the A/D converter 22 and encodes It using 
any of a number of different digital encoding methods such as, for example, 
differential phase shift l<eylng (DPSK), minimum shift keying (MSK), 
5 quaternary phase shift keying (QPSK), and quadrature amplitude 

modulation (QAM). IVIultiplier 26 receives the encoded data signal from 
encoder 24 and multiplies the signal with a radio frequency (RF) carrier 
signal and a pseudo noise (PN) sequence to create the transmit signal that 
is transmitted into the channel 16. It should be appreciated that multiple 
10 layers of multiplication may be performed in the multiplier 26 to achieve the 
transmit signal using the three Inputs. 

Multiplication by the RF carrier signal upconverts the encoded data 
signal to an RF frequency range centered about the earner frequency. DCO 
28 generates the RF carrier signal In response to a digital control signal from . 
1 5 the frequency hopping (FH) adaptation module 34 (and/or from the 
processing gain adaptation module 36). The N level chip sequence 
generator 30 and chip shaping filter 32 generate the pseudo noise (PN) 
sequence that is used to spread the transmit signal In frequency in a process 
known as direct sequence spread spectrum (DSSS). The N level chip 
20 sequence generator 30 produces a multi-level PN sequence in a manner 
that Is well known In the art. The chip shaping filter 32 filters the PN 
sequence In a manner that Insures low probability of intercept/low 
probability of detection (LPI/LPD). After filtering, the PN sequence exhibits a 
more Gaussian time domain distribution which makes hostile detection and 
25 parameter extraction more difficult. Further, the chip shaping filter 32 band 
limits the PN sequence to the channel, which reduces adjacent channel 
Interference. Note that a preferred embodiment of the invention can include 
frequency hopping (FH) adaptation and processing gain (PG) adaptation in 
concert, exclusively, or no adaptation. Without FH or PG adaptation, no 
30 feedback path Is required. 
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The FH adaptation module 34 is operative for implementing adaptive 
frequency hopping spread spectrum tecliniques in the system 10. The FH 
adaptation module 34 delivers a digital signal to DCO 28 instructing It to 
change the center frequency of the carrier signal based upon a frequency 
5 hopping sequence. By periodically changing center frequency, a number of 
advantages are obtained, such as increased resistance to jamming. In the 
preferred embodiment, the FH adaptation module 34 receives feedback from 
the receiver 18 which indicates, among other things, the type(s) and spectral 
location(s) of interference that currently exists in the channel 16. The FH 
to adaptation module 34 can detenmine from the feedback signals received 
from the Interference suppression processor 42 that a particular frequency 
hopping band is occupied by Interference that is difficult and/or impossible to 
overcome. The FH adaptation module 34 can then remove that frequency 
band from the frequency hopping sequence via frequency masking. When It 
5 Is detemiined that the frequency band is no longer in an unsatisfactory 

condition, the FH adaptation module 34 can again enter the frequency band 
into the frequency hopping sequence. In the preferred embodiment, both 
the FH adaptation module and the N level chip sequence generator 30 are 
synchronized by a sync signal 35 in a manner that is well known In the art. 
0 The processing gain adaptation module 36 Is operative for adaptively 

controlling the level of spreading that is Imparted to the modulated carrier 
signal. In this regard, the processing gain adaptation module 36 delivers 
control signals to both the N level chip sequence generator 30 and the chip 
shaping filter 32 to set the level of spreading. Like the FH adaptation 
5 module 34, the processing gain adaptation module 36 receives feedback 
signals from the receiver 18 that are indicative of the type(s) and spectral 
location(s) of interference presently in the channel 16. The processing gain 
adaptation module 36 uses this Infomiatlon to determine an appropriate 
processing gain for the transmit signal in light of the identified interference. 
3 For example, if a narrow band interference signal Is detected in the channel 
16, the processing gain adaptation module 36 may decide to spread the 
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modulated carrier an additional amount to further reduce the effect of the 
interfering signal. On the other hand, if little or no interference is detected in 
the channel 16, the processing gain adaptation module 36 may decide to 
reduce the level of processing gain (i.e., reduce spreading). 
5 In a prefenred embodiment of the present invention, the channel 16 is 

a wireless RF link. It should be appreciated, however, that the principles of 
the present invention can be implemented in communication systems having 
virtually any type of channel, including both wired and wireless channels. 
As illustrated in FIG. 2, the channel 16 includes an additive noise source 38 

10 and a source of interference/jamming signals 40. Additive noise source 38 
represents all sources of random noise in the channel 16. That Is, additive 
noise source 38 includes, for example, sources of themial noise, 
atmospheric noise, galactic noise, and others. Frequency selective fading or 
other types of fading can also occur in the channel 16. Interference source 

15 40 represents the sources of all man-made interference in the channel 16. 
For example, interference source 40 can include partial band noise 
jammers, spread spectrum co-channel sources, wideband noise jammers, 
bauded co-site signal sources, chirp jammers, multitone jammers, non- 
bauded co-site signals, and others. Signals from both additive noise source 

20 38 and interference source 40 combine with the transmit signal in the 
channel 16 to produce a modified signal. The modified signal is then 
received by the receiver 18. 

The receiver 18 includes: interference suppression processor 42, co- 
processor module 44, jammer signal of interest (JSOl) memory 46, 

25 conjugate and complex multiply (CCM) unit 48, N level chip sequence 
generator 50, chip shaping filter 52, fast Fourier transfomn (FFT) unit 54, 
inverse FFT unit 56, detector 58, decoder 60, post processor 62, digital to 
analog converter (DAC) 64, and performance estimator 66. As will be 
described in more detail, the interference suppression processor 42 

30 receives the modified signal from the channel 16 and suppresses the 
interference/jamming components therein using targeted interference 
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suppresslon. Co-processor module 44 provides supplemental processing 
power (such as, e.g., to perfom complex signal arithmetic) to the 
Interference suppression processor 42 for Increasing processing speed and 
thereby Improving data throughput. JSOl memory 46 stores, among other 
5 things, a library of feature plane representations corresponding to different 
jamming and interference types that may be encountered in the channel 16. 
Interference suppression processor 42 uses the information stored in the 
JSOl memory 46 to identify and classify interference components in the 
receive signal. As will be described In more detail, the Interference 
1 0 classifications are used by the Interference suppression processor 42 to 
detemriine an appropriate fonn of interference suppression for the receive 
signal. 

After the Interference/jamming signals In the receive signal have been 
suppressed, the interference suppression processor 42 outputs the restored 
15 spread spectrum signal for despreading, demodulation, and decoding. 
Con'elation of the output signal (i.e., to remove the pseudo noise 
modulation) is performed in the frequency domain (although other 
embodiments can use time domain con-elation). The N level chip sequence 
generator 50 and the chip shaping filter 52 are operative for creating the 

20 same PN sequence that was created In the transmitter 14 by N level chip 
sequence generator 30 and chip shaping filter 32. Operation of the N level 
chip sequence generator 50 Is synchronized using sync signal 100, which is 
a clock signal that is synchronized with the transmitter sync, signal, as will be 
familiar to those of skill In the art. The output of the chip shaping filter 52 is 

25 delivered to FFT unit 54 where it is converted to a frequency domain 

representation. Before conversion, the FFT unit 54 applies weighting to the 
pseudo noise sequence to reduce FFT edge effects. 

COM 48 receives the frequency domain output signal from the 
interference suppression processor 42 and the frequency domain PN 

30 infomnation from FFT unit 54. The COM 48 determines the complex 

conjugate of the data from the Interference suppression processor 42 and 
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then convolves this conjugate data with the frequency domain PN 
information. For reasons that are apparent to a person of ordinary sl<ill in the 
art, a significant time uncertainty can exist between the PN sequence 
produced in the N level chip sequence generator 50 and chip shaping filter 
5 52 and the PN sequence used to modulate the transmit signal in the 
transmitter 14. Use of frequency domain correlation in the receiver 18, 
rather than time domain correlation, provides significant computational 
savings when sample timing uncertainty is large. The convolved signal is 
delivered to inverse FFT 56 to be converted from a frequency domain 

1 0 representation to a time domain representation. 

Detector 58 receives the time domain correlated data from the inverse 
FFT 56 and detects significant correlation energy to produce detected data. 
In the preferred embodiment of the invention, detector 58 uses ensemble 
integration to perfomn the detection function. That is, detector 58 combines 

15 the magnitude squared of the inverse FFT 56 output. At zero or negative 
signal to noise ratio (SNR) conditions, symbols must be combined in an 
ensemble fashion and compared to a noise based threshold. The threshold 
is computed and updated dynamically during non-transmission periods in 
order to maintain an apriori bit error rate (BER) relative to ambient noise. 

20 Further processing may also be perfomied to adjust threshold levels during 
minimal background interference conditions. 

Following detection, the detected data is decoded and/or synthesized 
in decoder 60 to produce reconstructed data. Decoder 60 can perform 
adjustments to the multilevel code timing or phase in order to maximize the 

25 correlation peak of the signal. In this manner, the multi-symbol buffering 
nature of the detection process compliments blockwise, processor based 
frequency tuning offset correction methods, further improving correlation 
results. Decoder 60 produces in-phase and quadrature values from the 
detected data which are used to compute an instantaneous phase angle. 

30 The instantaneous phase angle is dealiased in the decoder 60 and 

decoding is performed via symbol based integration in order to produce the 
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reconstructed data. The reconstructed data is then passed to the post 
processor 62 where operations such as adaptive post filtering enhancement 
or frequency de-emphasis may be performed, producing conditioned, 
reconstructed data. The conditioned reconstructed data is delivered to the 
5 DAC 64 where it is converted to an analog signal for output to output device 
20. 

Performance estimator 66 receives the reconstructed data from the 
decoder 60 and analyzes the data to calculate one or more performance 
metrics. These performance metrics are then transferred back to the 
10 interference suppression processor 42 for use In fine tuning the interference 
suppression function. In a prefen-ed embodiment, signal to noise ratio 
(SNR), bit error rate (BER). and spectral distortion (SD) are used as 
performance metrics, although other metrics may also be used. 

As described above, the Interference suppression processor 42 is 

1 5 used to determine the type of interference that is present in the receive 

signal and to perfonm interference suppression on the signal based on the 
types of interference identified. The Interference suppression processor 42 
delivers signals back to the FH adaptation module 34 and the processing 
gain adaptation module 36 in the transmitter that are indicative of 

20 interference type and spectral location. Interference suppression processor 
42 Includes: frequency hopping tuner 68, FFT unit 70, signal memory 72, 
disturbance classification unit 74, jammer parameter extraction unit 76, and 
targeted interference suppression unit 78. Frequency hopping tuner 68 
receives the receive signal and processes the signal based on the known 

25 frequency hopping sequence of the transmitter 14. Frequency hopping 
tuner 68 also receives sync signal 100 to aid in signal processing. 
Frequency hopping tuner 68 outputs a time domain data signal that includes 
both desired information from the transmitter 14 and 
noise/interference/jamming signals from the channel 16. In an altemate 

30 embodiment, frequency hopping is not perfomned In the transmit unit 14 or 
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the receive unit 18 and the FH tuner 68 and the FH adaptation module 34 
are not required. 

The output signal from the frequency hopping tuner 68 is weighted (to 
minimize FFT edge effects) and processed In FFT unit 70 to convert It to a 
5 frequency domain representation. The frequency domain representation is 
then stored In the signal memory 72 for later use. Disturbance classification 
unit 74 receives both the frequency domain Information from FFT unit 70 and 
the time domain information from the frequency hopping tuner 68 and uses 
the infomrjation to identify and classify the interference/jamming components 

10 in the input receive signal. Classification Infonnatlon from the disturbance 
classification unit 74 is then delivered to the jammer parameter extraction 
unit 76 which determines parameters describing the Interference/jamming 
components based upon the classifications. That Is, the parameter 
extraction method (I.e., the types of parameters that are extracted and the 

1 5 method of calculating them) that Is used In the jammer parameter extraction 
unit 76 depends on the type of interference identified. By tailoring the 
parameter extraction method to the type of interference, the jammer 
parameter extraction unit 76 is able to detennine signal characteristics that 
are specific to the interference (e.g., baud rate, bandwidth, deviation, etc.). 

20 The targeted Interference suppression unit 78 Includes a library of 

software modules that are each capable of suppressing or removing 
undesired interference components from a subject signal. Each of the 
modules In the targeted interference suppression unit 78 wori<s best with a 
particular type or class of interference/jamming. For example, one module 

25 may be best at suppressing spread spectrum co-channel signals while 

another Is best at suppressing the effects of frequency selective fading. The 
targeted Interference suppression unit 78 receives the classification output 
signal from the disturbance classification unit 74 and the parameter signals 
from the jammer parameter extraction unit 76 and uses this Information to 

30 determine which module Is best for suppressing the interference in the 

receive signal. If more than one type of interference is present in the receive 
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signal, the targeted interference suppression unit 78 may choose more than 
one module to process the signal. After appropriate modules have been 
chosen by the targeted interference suppression unit 78, the unit 78 
retrieves the frequency domain and/or time domain receive signal 
5 information from the signal memory 72 and processes it using the 

appropriate modules. After processing, targeted interference suppression 
unit 78 outputs a frequency domain spread spectrum signal having 
suppressed interference. This signal is processed in the remainder of the 
receiver unit 18, as described above, to produce a data signal which is 

10 delivered to output device 20. 

As is apparent from the above description, the system 10 illustrated in 
FIG. 2 provides a hybrid system for interference mitigation. That Is, three 
distinct mitigation subsystems (targeted interference suppression, FH 
adaptation and processing gain adaptation) are implemented to handle 

15 virtually any type of Interference that can occur. The interference 

suppression processor 42 provides a receiver based pre-processor function 
for the removal of undesired jamming signals, effectively optimizing bit en-or 
rate perfomnance without the sacrifice of bandwidth or system capacity. The 
interference suppression processor 42 targets narrow band and partial band 

20 jammers such as partial band noise jammers, spread spectrum co-channel 
signals, bauded co-site and ambient signals, chirp jammers, multitone 
jammers, and non-bauded co-site and ambient signals. The FH adaptation 
module 34 provides a supplementary means of interference mitigation via 
interference masking and/or randomization of signal errors. Processing gain 

25 adaptation module 36 provides additional anti-jam protection by means of 
chip rate adjustment. Chip rate adjustments are especially effective in the 
presence of wideband jammer or co-channel spread spectrum interference 
that cannot be canceled or suppressed via other methods. In applications 
where size, weight, power, and bandwidth limitations are important, such as 

30 handheld or manpack systems, the FH adaptation and processing gain 

adaptation components of the system 10 can be omitted. Such an omission 
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removes the need for the FH adaptation module 34, the processing gain 
adaptation module 36, the feedback channel from the interference 
suppression processor 42, the frequency hopping tuner 68, and the digitally 
controlled oscillator 28 (which can be replaced by a single frequency 
oscillator). 

The disturbance classification unit 74 processes the time and/or 
frequency domain receive signal data to produce a maximum likelihood 
estimate of the disturbance class superset and/or subset classes. Class 
supersets can include narrow band, partial band and wideband classes. 
Class subsets, on the other hand, can include single tone, multi-tone, swept 
tone, bauded communication signal, non-bauded communication signal, 
partial band noise, or wideband noise. Interference may also be classified 
by modulation type, such as, for example, FSK, NRZ-ASK, and BPSK 
modulation. The disturbance classification unit 74 can utilize any of a 
number of techniques to recognize jammer and/or interference signals. For 
example, the disturbance classification unit 74 can employ tree based 
classification, multi-layer perceptron classification (MLP), or other methods 
of classification. 

FIG. 3 is a block diagram illustrating the disturbance classification unit 
74 in one embodiment of the present invention. As illustrated, the 
disturbance classification unit 74 includes a feature extraction unit 80 and a 
classifier 92. The feature extraction unit 80 receives time domain and/or 
frequency domain representations of the receive signal and processes this 
data to derive certain features and parameters that are unique to the 
interference/jamming components within the receive signal. These features 
and parameters are then delivered to the classifier 92 for use in classifying 
the interference and jamming components. In the illustrated embodiment, 
the feature extraction unit 80 includes: super class recognition unit 82, 
adaptive super class masking unit 84, subclass feature extraction unit 86, 
subclass interference suppression unit 88, and subclass parameter 
extraction unit 90. Superclass recognition unit 82 is operative for deriving 
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parameters that are useful in classifying the interference/jamming 
components according to superclass. The superclass masking unit 84 is 
operative for signal separation, for example, to separate FSK, ASK and PSK 
signals. The subclass feature extraction unit 86 performs feature plane 
5 transfomiations on the receive signal to aid in the classification of 
interference^amming components based on subclass. Subclass 
interference suppression unit 88 Is operative for SNR improvement for the 
isolated signal. Subclass parameter extraction unit 90 determines subclass 
parameters from the input data based on the features identified by the 
10 subclass feature extraction unit 86. 

The classifier 92 receives the features and parameters determined in 
the feature extraction unit 80 and uses them to classify the 
interference/jamming components in the receive signal. To expedite the 
classification process and facilitate real-time operation, a separate co- 
15 processor module 44 can be utilized by the classifier 92 and block 80. As 
part of the classification process, the classifier 92 compares features and 
parameters detennined by the feature extraction unit 80 to those of known 
Interference/jamming types that are stored in the JSOl memory 46. Use of a 
J SOI memory 46 allows new interference types to be added to the system in 
20 a quick and efficient manner in contrast to prior art methods that are fixed, 
"hard wired" implementations. In one embodiment of the invention, a 
content addressable memory (CAM) is utilized to both store the known 
interference types and to perform the required comparisons. As is well 
known in the art, a CAM is capable of comparing an input signal to the 
25 contents of a plurality of memory locations simultaneously. The CAM then 
outputs a signal indicating all memory locations that include the input signal. 
The classifier 92 outputs a signal that is indicative of the type of interference 
in the receive signal. Such classification can include, but is not limited to, 
the following bauded and unbauded signals of interest: 
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Amplitude Modulation 
PAM (Pulse Amplitude Modulation) 
ASK (Amplitude Shift Keying) 

5 Phase Modulation 

M-ary PSK (Phase Shift Keying) 

Offset M-ary PSK 
PPM (Pulse Position Modulation) 

"•0 Frequency Modulation 

M-ary FSK (Frequency Shift Keying) 
MSK (Minimum Shift Keying) 
FDM/FM (Frequency Division Multiplexing/ 
Frequency Modulation) (non-bauded) 

Amplitude/Phase Modulation fAM/PM> 
M-QAM (Quadrature Amplitude Modulation) 

In one embodiment of the present invention, the superclass and subclass 
20 feature extraction and interference suppression functions are combined into 
a single unit and all recognition functions are performed by the disturbance 
classification unit 74. 

In the preferred embodiment of the invention, the feature extraction 
unit 80 perfonms an efficient data analysis using circa 50 data transitions. 
25 For extremely low baud rates, feature extraction unit 80 also performs a 
feature accumulation function whereby features of contiguous signal 
sources are saved until sufficient data transitions have been accumulated for 
identification. In a preferred embodiment of the invention, a general- 
purpose set of classification parameters are computed by feature extraction 
30 unit 80. Classification parameters are computed in the feature extraction 
unit 80 using at least one of a plurality of feature plane transformations. In 



BNSDOCID: <VVO 9939445A1 I > 



-18- 



10 



15 



one embodiment, the feature plane transfomnations are computed from a 
joint time-frequency (JTF) matrix H of order n.m wliere n represents a 
contiguous time index and m represents a contiguous spectral index. 

FIG. 4 illustrates a JTF feature plane showing three discrete 
interference components versus time. The first component 102 corresponds 
to an FSK signal, the second component 103 corresponds to an non retum 
to zero (NRZ)-ASK signal, and the third component 104 corresponds to a 
PSK signal. JTF matrix H quantifies the information in the JTF graph. The 
JTF matrix can be represented as follows: 



H - 



KO,i) 



h(0,n- 1) 



KIO) h(2,0). h(m- 1,0) 



h(m-l,n-l) 



Eq. 1 



The feature plane transformations are perfomied by using the JTF matrix H 
to compute a plurality of transfomiation matrices. The following are some 
transformation matrices, TN, that may be used in this regard (although other 
similar transformations can also be used): 



Transfomriatlon MatriV 1 



20 



T 1 =: 



a(0,l) fl(0;2) a(0,3) 
0 a(l,2) a(l,3) 
0 0 a(2,3) 



a(0,/w-l) 
a(l,w-l) 
a(2,m-l) 

a(fn-2,m-Y) 



Eq. 2 



where 
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n-l 



n-1 



a(U) 



^ h{l. i ) h (J, i) - ( X ( X ^" 

i=0 



w-1 



[ ( X h(J.i)^-(X NJ.i))^/n )( X *^(».*)'-(X h(U))Vn ) } 

1=0 /sO 1=0 



/I- 1 

I 

/=0 



Eq. 3 



T 2 = 



Transformation Matrix 2 

■Z»(0,1) 5(0,2) !>(03) 6(0.«-l) ■ 
0 6(1,2) 6(1,3) 6(1, n-l) 
0 0 6(2,3) 6(2, «-l) 

. 0 . ft(«-2,«-l) 



Eq. 4 



where 



m-l 



X h{U)h(U).(X I))(X J))/m 



/sO 



i=0 



15 



[ ( X h(i.J)^-(X h(i,J))2/m )( X h{U)*-(X Ni.l))*/m ) }« 

<=0 i=0 reO ieO 



Eq. 5 



20 



Transformation Matrix 3 



T3 = [p^, Pa, Pa ..p«.i]. Where p^ = max h(k,l)' V i 



Eq. 6 



Transformation Matrix 4 



25 T4 = [Ho, m J, here ^t^ = mean h{k,\y V i Eq. 7 
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^(0,0) 

dm) 



d(0,n - 1) 
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Transformation Matrix 5 
rf(l,0) </(2,0) dim -1,0) 



Eq. 8 



10 



where d(l,J) = 
d(l,J) = 



+1,h(I.J)2>c*P„ 
0, otherwise 



Transformation Matrix 6 



T6 = [Ho, \L^, Where = mean d(k,n) V n Eq. 9 



Transformation Matrix 7 



15 



20 



In addition to the transformation matrices set out above, 
computationally efficient signal mapping methods are used In accordance 
with the present invention to generate unique stationary patterns In order to 
supplement classification, while eliminating carrier frequency offset errors. 
Since the exact earner frequency is undetemiined before the mapping 
procedure, a balanced multiplication Is used in order to cancel the carrier 
components. Representing the input in-phase and quadrature signal 
samples as x(i) and y(i), and the transformation in-phase and quadrature 
samples as px(i) and py(i), the transformation is as follows: 



25 



pxG) = y(i+t)[y(l)-x(i)] + x(i+T)[x(i) + y(i)] Eq. 10 
PVO) = x(i+T)[x(i)-y(i)] + y(i-H:)[y(i) + x(i)] Eq. 1 1 



where: 
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T = delay, in samples 
i = 1 to N 
j = 1 toN-x 

5 For PSK signal formats, the desired delay t is equal to the bit period. 

By up sampling to achieve a 10 to 1 oversamplrng, a close approximation to 
the desired delay can be obtained. Projection of the histogrammed polar 
space defined by the above transformation to a Cartesian space S gives the 
following transformation matrix which defines the unique signal signature, 
10 where L and Q are arbitrary integer constants corresponding to magnitude 
(L) and phase (Q) dimensions. 



T7 = 



5(0,0) 
S(0,l) 
S(0,2) 



5(1.0) 5(2,0) SiQ-1,0) 



SiQ-l,L-l) 



Eq. 12 



15 T7 projections provide additional feature transformation data, as 

shown in Equations 13 and 14 below. 



Transformation Matrix 8 



20 T8 = ^i,, ^2 J, where = mean S(k,l) V I Eq. 13 

Transformation Matrix 9 

T9 = [^o> m n,^.,], where ii^ = mean S(q,k) V q Eq. 14 

25 

For FSK signal formats, the desired delay t is not necessarily equal to 
the bit period. The pattern is not stationary as in the PSK case where the 
pattern peaks remain relatively fixed in the phase coordinate dimension. 
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The changing characteristic of the pattern structure for FSK formats can be 
used to distinguish signal fomiats that have phase slopes versus fixed 
phase states during the bit periods. This characteristic is made part of the 
signal descriptive vector and used to identify FSK formats. For pattern 
structures that change, phase derivative processing Is used to derive a 
second signal pattern. 

In order to obtain the optimum delay for a PSK signal, a baud rate 
estimate Is employed. A delay of approximately one half bit, as estimated by 
the pre-filter process spectrum width detennination. Is used. The actual 
baud rate is determined by perfomiing a spectnim analysis of the transform 
output. From this baud rate, the proper one bit delay time can be 
determined. Methods for baud rate calculation are discussed below. 

FIG. 5 illustrates a cross-section of a typical feature plane signal 
transformation for a five FSK interference signal, obtained using equations 
10, 1 1, 12, and 13 above. It should be appreciated that the full feature plane 
transformation is a three dimensional representation and the two 
dimensional cross-section of FIG. 5 is presented for ease of illustration. 
Feature plane transfonnatlons, such as the one In FIG. 5, de-emphasIze the 
random components of the signal and emphasize the invariant features, 
supporting accurate classification of signal type by producing a unique 
signal signature. Note that the transformation includes five distinct frequency 
modes that aid In identification of Interference type. A plurality of feature 
plane transformations representing different interference types can be stored 
In the JSOl memory 46 for comparison with transformations derived from 
receive signals to classify interference components in the receive signals. 

Since the signal processing required to perform the signal 
transformations includes non-linear processing, It is desirable to implement 
as much signal to noise ratio (SNR) Improvement ahead of the signal 
processing as possible. In one embodiment of the Invention, SNR 
Improvement is perfomied before non-linear processing in the H transform 
domain, as well as after non-linear processing in the TN transform domain. 
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It is highly desirable that the SNR improvement not distort the signal features 
that produce the pattern for signal recognition. Furthemiore, It Is desirable to 
separate co-exIstent signals in order to minimize feature aliasing. Two 
methods of SNR processing that can be employed are preciassification 
5 filtering and phase domain filtering, in the prefenred embodiment of the 
present invention, both methods are used. 

FIG. 6 is a flowchart illustrating a superclass classification and 
masking method to perform preciassification filtering in accordance with one 
embodiment of the present invention. Prior to the operations performed by 

1 0 FIG. 6, it is assumed that the H matrix and transformations T1 through T9 
have been computed. Superclass classification begins with Select Min T1 
260, which searches the anray to produce the smallest correlation. Next, 
ortiiogonal 262 decision compares the value produced by Select Min T1 to 
an apriori negative constant to determine If a significant negative correlation 

15 Is present. If so, decision block 264 examines T5 for detection. If detection 
Is positive nonzero, Mask SCO 266 then determines signal bandwidth from 
T1 indices n and m, which represent the signal extent for Super Class O. As 
a further step, class SCO may be confirmed by comparing the peak power 
from columns n and m for consistency. Finally, Mask SCO 266 filters the 

20 SCO signal by zeroing the identified band, and stores the mask location to 
memory. Following consideration of all T1 (i.e., ail SCO), the Superclass 
Classification continues with Select T6 SS 270, which examines the T6 
transfomnation for near steady state (SS) values by comparison against an 
apriori constant. For SS T6, the maximum signal energy is determined from 

25 Find Max T4, which coarsely locates the center frequency of the candidate 
signal, producing a peak index. Next, decision block 274 examines T5 for 
detection. If detection is positive nonzero, Orthogonal Proximity 276 
compares T1 (corresponding to the index from Find Max T4 and proximity 
indices) against a second apriori negative constant to determine if significant 

30 negative correlation is present. If so, then Mask SCI determines the 
bandwidth of the signal bounded by the contiguous detected SS values 
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about the max T4 index, and filters the SC1 signal by zeroing the identified 
band, and stores the SC1 mask location to memory. If significant negative 
correlation is not present, then Mask SC2 280 determines the bandwidth of 
the signal and filters the SC2 signal by zeroing the Identified band, and 
5 stores the SC2 mask location to memory. Next, if ail distinct steady state 
regions have been considered (i.e., all candidate signals within the analysis 
band), then ail T6 SS 282 branches to End. In this discussion SCO 
represents an FSK signal superclass, SCI refers to bauded signal 
superclasses (e.g., PSK), and SC2 represents non-bauded continuous 

10 wave (CW) signals. The masking and classification process of FIG. 6 serves 
as preliminary signal separation and recognition prior to sub-class 
parameter extraction and final classification. 

After signal separation and superclass classification, the following 
sub-class parameters may be computed, although others may also be 

15 appropriate: 



1 


Dynamic Range, 


DR 


2 


Amplitude Modes, 


AM 


3 


Amplitudes RMS, 


AR 


4 


Discriminator Deviation, 


DD 


5 


Phase Modes, 


PM 


6 


Discriminator Mean, 


DM 


7. 


Outlier to Mean Ratio. 


OM 


8 


Outlier to RMS Ratio, 


OR 


9 


Baud Rate, 


BR 


10 


Superclass, 


SC 


11 


Frequency Bandwidth, 


BW 


12 


Frequency Modes, 


FM 


13 


Frequency Center, 


FC 



30 
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Transformations T7, T8, and T9 may be used to compute mode information 
by tabulation of peak infonnatlon, for example, using derivatives of filtered 
transformations. Alternatively, the modes may be determined in a similar 
fashion using filtered histograms of dealiased phase Infomriation, magnitude 
5 information, or discriminator output. Discriminator output also provides 
signal bandwidths and center frequencies. The above-described features 
provide sufficient information for classification into multiple signal types 
including, but not limited to: CW, FSK, MPSK, QPSK, MASK, and FM 
Unbauded. It should be appreciated that preclassificatlon filtering can be 

10 used to perform functions other than signal separation. 

In one embodiment of the present invention, preclassificatlon filtering 
is performed by isolating the narrow band interference signal of interest, 
estimating the bandwidth of the signal of interest, and prefiltering the signal 
of interest before classification functions are invoked. FIGS. 7-1 1 illustrate 

15 this method. First, as illustrated in FIG. 7, the signal of interest is converted 
into a frequency domain representation using a frequency transformation, 
such as a discrete Fourier transform (DFT) or a fast Fourier transform (FFT). 
The spectrum illustrated in FIG. 7 represents an MSK signal. The frequency 
coefficients found using the frequency transformation are converted to power 

20 spectrum and smoothed by a moving average filter, or similar device, to 
achieve the smooth spectrum illustrated in FIG. 8. The smooth spectrum is 
histogrammed, as illustrated in FIG. 9, to determine the noise floor level of 
the interference signal of interest (i.e., the peak of the histogram). Using the 
peak spectrum location from the histogram as the location of the signal noise 

25 floor 96, the signal bandwidth at the noise power level is measured as 

illustrated in FIG. 10 producing bandwidth 98. This width is corrected to a 
null-to-null bandwidth by assuming the signal spectrum shape as sin(x)/x 
and utilizing the peak SNR reference to determine where on the sin{x)/x 
curve the width was measured. FIG. 1 1 illustrates the spectrum of the MSK 

30 signal after the prefiltering operation. 
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As described previously, phase domain filtering is another technique 
that can be used to improve SNR before nonlinear processing. Phase 
domain filtering involves filtering the de-aliased phase components of the 
signal of interest using a moving average filter although other filter 
embodiments can also be used. The moving average filter that is used 
should possess sufficient bandwidth to prevent distortion of the signals of 
Interest. In one embodiment of the present invention, a three point moving 
average filter was found to possess sufficient bandwidth. The following 
equation describes one embodiment of a discrete time averaging filter that 
may be used to perfomi phase domain filtering: 

The moving average filter is preferably applied repeatedly to the signal of 
Interest. Phase domain filtering has been found to be effective on FSK. 
MSK, and PSK signals. Use with other types of signals is also possible. 

FIGS. 12 and 13 are cross-sections of feature plane signal 
transfomiations illustrating the beneficial effects of using phase domain 
filtering. FIG. 12 Illustrates the T7 transform projection that is obtained using 
a five tone FSK signal at 1 dB SNR (after spectral excision interference 
suppression, as will be described shortly). Such a transform should display 
five distinct feature peal<s. As illustrated in FIG. 12, these peaks are not 
present with enough signal power to discern the presence of the five tone 
FSK signal. FIG. 13 illustrates the feature plane transform for the same 
signal after phase domain filtering is applied. Note that the five distinct 
frequency peaks are clearly visible. As can be appreciated, use of the 
transformation illustrated in FIG. 13 significantly enhances the likelihood of 
proper signal classification. The phase domain filtering has effectively 
increased the SNR of the transform infomriation by approximately 6 to 7 dB. 
Classification and baud rate procedures can then be used to process the 
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improved transform. In order to overcome the noise enhancement effects 
introduced by derivative computation, phase domain filtering may be 
perfomned prior to application of the discriminator. 

As described previously, the baud rate of the interference signal of 
5 interest can be one of the extracted features that is delivered to the classifier 
92 for use in classifying thie signal of interest, in accordance with the present 
invention, methods are provided for determining the signal baud rate of the 
interfering signal of interest based on the modulation method of the signal of 
interest. That is, first it is necessary to determine the type of modulation 

10 being used (e.g., frequency modulation or phase modulation) and then a 
baud rate detemilnation method is chosen based on modulation type. If an 
interfering signal of interest is found to have phase modulation, baud rates 
are detemnined by analyzing the frequency content of the T7 feature space 
of the signal of interest represented by equations 10 and 1 1 . That is, first a 

15 T7 feature space transformation is found for the signal of interest and then 
the feature transformation is converted to a frequency domain representation 
using, for example, an FFT. The frequency domain representation is then 
analyzed to isolate a rate line corresponding to the baud rate. Even without 
prefiltering, a rate line can be detected at low SNRs (i.e., Odb and below) 

20 using this method. FIG. 14 is a graph illustrating the spectrum of the T7 
feature space of a BPSK signal at OdB SNR. Note that the rate line 94, 
corresponding to 250 kHz, is clearly visible in FIG. 14 and may be readily 
extracted by a simple search algorithm. 

For frequency modulation (FM) signals, a rate line may be established 

25 by obsen/ing the transitions in the T7 space. The transitions are used to 
generate a linear array of alternating Kronecker deltas as illustrated in FIG. 
15 in connection with a five tone FSK signal. A frequency transfonnation of 
the transition sequence produces a frequency spectrum (see FIG. 16) that 
clearly indicates the desired baud rate line 95. 

30 As described previously, the targeted interference suppression unit 

78 includes a plurality of interference suppression modules that can be used 
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to suppress interference/jamming components within the receive signal. In 
accordance with the present invention, a plurality of interference 
suppression techniques are provided that can be implemented as modules 
within the targeted interference suppression unit 78. These techniques are 
(1) the inverse whitening function technique, (2) the adaptive Inverse weight 
technique, and (3) the adaptive excision technique. In addition to these 
Interference suppression techniques, as described previously, the system 10 
also provides frequency hopping adaptation and processing gain adaptation 
functionality to help it adapt to hostile spectral environments. The above 
techniques will now be described. 

In a preferred embodiment of the invention, first and second order 
statistics of the receive signal spectrum are derived prior to application of 
Interference mitigation methods. That Is, non-arlthmetic modal moment 
estimates are used to describe different "modes" In the spectrum. Each 
mode corresponds to a single component in the receive signal including the 
fundemental communication signal and interference components. Modal 
moment estimates avoid jammer induced power bias that can arise when 
using arithmetic estimates. Modal moment estimates also allow multlmode 
spectra to be accurately described and ultimately suppressed. Modal mode 
estimates are also tolerant of wide dynamic range fluctuations. 

FIG. 17 Is a plot Illustrating typical modal moment estimates for a data 
signal corrupted by a partial band noise jammer. The x axis of the plot 
represents the spectral magnitude of the receive signal. The y axis of the 
plot represents the frequency with which the receive signal achieved a 
particular magnitude. That is, the y axis Is related to the number of samples 
In the Input data that are within a particular magnitude range. As such, the 
plot represents a magnitude distribution for the frequency domain 
representation of the signal. As shown in FIG. 17, two modes are apparent 
from the plot. The first mode corresponds to the fundamental data signal 
(l.e., the desired infomiation from the transmitter). The second mode 
con-esponds to the partial band noise jammer. In accordance with one 
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embodiment of the present invention, these two modes are described using 
modal estimates p. and ct, wherein }i is the mean magnitude and o is the 
magnitude standard deviation of the mode. An adaptive limit a can also be 
determined from the plot to describe a limitation point for the first mode. In a 
5 preferred embodiment of the present invention, the adaptive limit a for the 
fundamental data signal is equal to ka+|i, where k is a constant. 

FIG. 18 is a flowchart illustrating a process for deriving the modal 
moment estimates from the receive signal. A signal is first received from the 
channel (step 130). The receive signal is transformed to a frequency 

10 domain representation using a frequency transformation, such as an FFT 
(step 132). Magnitudes of the frequency transformed signal are then 
computed (step 134) and a modal envelope is developed (step 136). The 
modal envelope is then filtered (step 138) and a derivative is taken to aid in 
the detection of modes within the envelope (step 140). For example, the 

1 5 derivative will equal zero at the peaks of the modal envelope. Once the 
modes are identified, the mean magnitude |x of each mode is located (step 
142). A standard deviation c is then measured for each mode (step 144). 
When modal moment estimates for each of the modes have been computed, 
an adaptive limit a is calculated (steps 146 and 148). The modal estimates 

20 are then stored for use in interference suppression (step 150). 

As described above, in one embodiment of the present invention the 
adaptive limit a is calculated as the sum of (1) the mean magnitude of the 
fundamental data signal n,, and (2) the product of a constant k and the 
magnitude standard deviation a, of the fundamental data signal. In one 

25 embodiment of the invention, the constant k can be changed based on 

various factors, such as the type of interference that is present. For example, 
the constant k can be changed based on an estimate of relative jammer 
signal strength and bandwidth in order to obtain near optimum bit error rate 
performance. 
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Before discussing interference suppression metliods. an expanded 
descriptfon of FH and processing gain adaptation will be made. FIG. 19 is a 
spectrum diagram illustrating a number of bands WsN for use in frequency 
hopping. Although only four bands are shown, it should be appreciated that 
N can be virtually any number depending on the application. The bands 
WsN ail occur within a communication system bandwidth Wt. As illustrated, 
each band WsN includes an Interference signal. That is, band Ws1 includes 
a frequency selective fade 110, band Ws2 includes a partial band noise 
jammer 112, band Ws3 includes multiple tone jammers 114, and band Ws4 
includes bauded Interference 116. As described above, the frequency hop 
adaptation module 34 will sequence through a predetemnined hop 
sequence during normal operation, wherein the frequency of operation of 
the system 10 will alternate among bands WsN. In accordance with one 
aspect of the present invention, a particular band can be removed from the 
frequency hop rotation based on the interference content within the band. 
For example, the system 10 can determine that effective communications 
cannot be carried out in band Ws2, due to the presence of the relatively 
powerful noise jammer 1 12, and Ws2 is removed from the frequency 
hopping sequence. 

FIG. 20 is a flowchart illustrating one method for perfomriing frequency 
hop masking. First, the transmitter 14 hops to the next frequency band in the 
hop sequence (step 152). The transmitter 14 then transmits a signal into the 
channel 16 within the new frequency band (step 154). For example, the 
transmitter 14 may transmit a spread spectrum signal within band Ws2 
illustrated In FIG. 19. The signal is received from the channel 16 by the 
receiver 18 and is processed to detemriine interference classification (step 
156). Modal estimates \i and a are then determined for each of the modes 
identified (step 158). The modal estimates are then used to detennine 
estimates for interference power and interference bandwidth for each of the 
identified modes (steps 160 and 162). Next, the signal is classified (step 
163). An estimate of the bit error rate that will be achieved by demodulating 



I > 



wo 99/39445 



PCT/US99/01897 



-31- 

and decoding the receive signal without interference suppression Is 
calculated using the estimated power and bandwidth infomiatlon and the 
classification infornnation (step 164). The estimated bit error rate is then 
compared to a maximum tolerable bit error rate BER^^^ (step 166). If the 
5 BER does not exceed the maximum tolerable bit error rate, the system 1 0 
hops to the next frequency band and the process is repeated. If the BER 
exceeds the maximum tolerable BER, the present hop band WsN is masked 
from the hop sequence so that it is no longer used by the frequency hopping 
adaptation module 34 (until the interference condition changes) (step 170) 

1 0 and the system 1 0 then hops to the next frequency band. 

The procedure outlined In the flowchart of FIG. 20 can also be used 
for processing gain adaptation. That is, Instead of changing the hop 
sequence in step 170, the processing gain parameters determined by the 
processing gain adaptation module 36 can be changed based on the 

15 identified interference signals. In general, new processing gain parameters 
will be chosen to minimize the effect of the Identified interference. The 
decision on whether to use frequency hopping adaptation or processing 
gain adaptation will depend on both the type or classification, spectral 
location, and level of the interference identified by the disturbance 

20 classification unit 74 and the jammer parameter extraction unit 76 and the . 
desired performance goals of the system 10. Circumstances may exist, for 
example, when both frequency hopping adaptation and the processing gain 
adaptation are performed. 

In a preferred embodiment, the above-described FH/processIng gain 

25 adaptation procedure Is performed every time a receive signal is received 
from the channel 16. Alternatively, as will be described shortly, the 
procedure will only be performed if the selected interference suppression 
method does not achieve desirable results. In such an embodiment, the 
interference modules will access an FH/processing gain adaptation module 

30 if a predetermined perfomiance goal (such as a predetemriined BER) is not 
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achieved. A description of the Interference suppression methods of the 
present invention is now made. 

INVERSE WHITENING FUNCTION 

5 

The inverse whitening function suppression method makes use of the 
fact that direct sequence spread spectrum data transmissions generally 
have a spectrum similar to that of band limited white noise. Based on this 
fact, perturbations in the channel serve to color or distort the spectral 

1 0 characteristic of the channel from this white noise model. The inverse 

whitening function removes this distortion from the white noise model. FIG. 
21 is a flowchart illustrating the operation of the inverse whitening function in 
one embodiment of the present invention. A signal is first received from the 
channel 16 (step 174). An estimate of the spectral envelope of all interfering 

15 signals within the receive signal Is then obtained (step 176). In order to 
insure that the underlying phase Information is undisturbed, an all-pole 
spectral envelope estimate is required. In one embodiment of the present 
invention, the following all-pole spectral model is implemented: 

Where a^ represents the prediction coefficients for an nth order model. The 
prediction coefficients may exemplify the periodic state estimate that is 
obtained by an adaptation processor using, for example, classic auto- 
25 correlation and durbin recursion. 

After the spectral envelope is estimated, the receive signal is inverse 
filtered using the inverse of the spectral envelope to remove the interference 
components (step 178). In one embodiment of the invention, the inverse 
whitening filter is represented by the following equation: 
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y{n) = x{n)+Y.a{k)x{n~k) Eq.16 

After whitening, the whitened spectrum is delivered to the CCM 48 to Initiate 
5 despreading and demodulation as described above (step 180). The Inverse 
whitening procedure is then repeated for the next receive signal. 

If the signal having the whitened spectrum still results in a relatively 
high BER, the system may want to perform FH or processing gain 
adaptation. In this regard, the following steps are also performed. Using the 

10 time domain residual signal and the receive signal, interference power 

estimates are calculated (step 184). The power estimates and interference 
classifiction are used to estimate the BER that will be achieved (step 186). If 
the estimated BER exceeds a maximum BER, then frequency hop adaptation 
or processing gain adaptation can be performed as described above (step^ 

15 188). BER estimates can also be computed using the modal approach of Fig. 
20. 

In an alternate emk>odlment, the inverse whitening approach is 
applied In the frequency domain. The Interfering signal is located using 
modal estimates and subsequenct noise floor characterization. The location 

20 of the interfering signal, the bandwidth of the signal, and the classification of 
the signal are used to apply an inverse filter envelope obtained from an 
apriori family of such envelopes. FIG. 22 is a spectmm diagram illustrating 
the spectrum of a spread spectrum data signal 120 corrupted by a PSK 
partial band jammer 122, In accordance with the invention, estimates are 

25 made of the spectral envelope of the PSK jammer 122 and an inverse 

matched filter estimate is derived (not shown). As Illustrated in FIG. 22, the 
bandwidth (BW) and the center frequency (fo) of the jammer 122 are 
estimated, and an Inverse sin(x)/x matched filter response (not shown) is 
generated. The receive signal is applied to the inverse matched filter to 

30 obtain the whitened output spectrum with no reduction in system throughput. 
The inverse whitening function Interference suppression method can be 
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used to suppress a wide range of narrowband and partialband interference 
types, including tones, communication signals, and jammers. In addition, 
the inverse whitening function method is capable of improving bit error rate 
performance at high interference power level to signal power level ratios. 
5 The inverse whitening function method can also be used to remove the 
damaging effects of frequency selective fading from the received signal. 

ADAPTIVE INVERRg WFIfiHT 

"•0 In the adaptive inverse weight interference suppression method, the 

modal moment estimates are used to selectively target all spectral outliers 
that exceed the adaptive modal limit a. As described earlier, the adaptive 
modal limit can be calculated as n+ka. Only those signals that exceed this 
limit are affected by the anti-jam suppression. For example, with reference 

15 to FIG. 17, only those points to the right of adaptive limit a on the graph are 
affected by the interference suppression method. An inverse weight © is 
applied to each of these samples that, in one embodiment, is inversely 
proportional to the observed deviation from a (I.e., inversely proportional to 
mag(sample) - a). In one embodiment of the Invention, the weights o) are 

20 determined by the following equation: 

fl) = a/[a-h(5-a)"'] Eq.17 

where m is an integer. 

25 The adaptive inverse weight method can also be applied to proximity 

samples surrounding each spectral outlier sample, thereby further improving 
performance by enhancing the attenuation of the Interference envelope. 
That is, an inverse weighting can be applied to +/- n samples on either side 
of each sample that exceeds the adaptive limit a. The adaptive inverse 

30 weight interference suppression method Is effective for reducing the BER for 
signals containing a myriad of different interference types, including 



BNSDOCIDkWO g93944SA1 I > 



wo 99/39445 



PCTAJS99/01897 



-35- 

intentionai hostile jammers, and bauded and non-bauded interference. 
Also, the adaptive inverse weight method is effective at very high jammer 
power to signal power ratios. 

FIG. 23 is a flowchart illustrating the operation of the adaptive inverse 
weight Interference suppression method. First, a signal is received from the 
channel 16 in the receiver 18 (step 194). The signal is then frequency 
transformed to achieve a frequency domain representation and modal 
estimates ji.a are computed for each of the identified modes (step 196). 
Next, an adaptive limit a is calculated using the modal estimates of the 
fundamental data signal (step 198). Each of the samples from the frequency 
domain infomnation Is then compared to the adaptive limit a (step 200). If 
the magnitude of a spectrum sample is greater than the adaptive limit, an 
inverse weight filter is applied to that sample (and optionally to proximity 
samples) as described above (steps 202 and 204). When all of the samples 
have been processed, the filtered signal is delivered to the CCM 48 to 
initiate despreading and demodulation (steps 206 and 208). Processing is 
then initiated for the next received signal. 

To determine whether to perform FH or processing gain adaptation, 
interference power is estimated using the received signal and the adaptive 
limited data (step 212). The BER is estimated (step 214) and, if the BER is 
greater than a maximum value, frequency hopping or processing gain 
adaptation is performed (step 216). 

Using the adaptive inverse weight interference suppression method, 
improved signal fidelity is obtained without sacrificing throughput. In general, 
the interference is suppressed without significant distortion of the underlying 
spread spectrum data. 

ADAPTIVE FXniRinN 

In the adaptive excision method, the modal moment estimates are 
used to locate and characterize the interference signal components. An 
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estimate of jammer signal bandwidth is tlien used to excise (i.e., remove) a 
portion of the signal spectrum from the receive signal, including both 
interferer and signal. Unlike the adaptive inverse weight method, the 
adaptive excision method removes a portion of the spread spectrum signal 
spectrum from below the adaptive limit a. Consequently, the adaptive 
excision method does not achieve the bit en-or rate reductions that the 
adaptive inverse weight method is capable of. However, the adaptive 
excision method is simpler to Implement than the adaptive inverse weight 
method. 

In one embodiment, a filter response modeling an ideal notch filter is 
used to remove the portion of the spectrum occupied by the interfering signal 
or signals. Other filter configurations, such as sinusoidal rolloff functions, 
can alternatively be used to minimize adverse affects during reverse 
frequency transformation caused by sharp spectral edges. As with the other 
methods, the adaptive excision method improves signal fidelity without 
sacrificing bandwidth, 

FIG. 24 is a flowchart illustrating one embodiment of the adaptive 
excision Interference suppression method. A signal is first received from the 
channel 16 (step 220) and the signal is processed to compute the modal 
estimates n, o (step 222). The adaptive limit a is then calculated using the 
modal estimates for the fundamental data signal (step 224). The frequency 
domain signal samples are then each compared to the adaptive limit a (step 
225). If the magnitude of a sample is greater than a, a notch like filter 
response is applied to the sample (steps 226 and 228). Once all of the 
frequency domain samples have been processed, the frequency domain 
Information is delivered to the COM 48 to Initiate despreading and 
demodulation of the data (step 230 and 232). As with the other methods, 
processing is also performed to determine whether frequency hop or 
processing gain adaptation is to be performed (steps 234-240). 
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As described previously, each of the above described interference 
suppression methods can be implemented in the targeted interference 
suppression unit 78. In addition, other interference suppression methods, 
such as prior art methods, can be implemented in the targeted interference 
5 suppression unit 78. The unit 78 will determine which rnethod to use based 
upon the type or types of interference detected In the disturbance 
classification unit 74 and the jammer parameter extraction unit 76. It should 
be appreciated that It is also possible to cascade suppression methods in 
accordance with the present invention. That is, two or more suppression 

1 0 methods can be implemented in the targeted interference suppression unit 
78 if it is determined that this is the best way to mitigate interference. FIG. 25 
is a flowchart illustrating one such cascaded interference suppression 
procedure. The signal Is received from the channel 16 in the receiver 18 
(step 240). An adaptive inverse weight method is applied to the receive 

15 signal to remove one type of detected interference from the signal (step 
242). An inverse whitening function is then applied to the receive signal to 
remove another type of Interference from the signal (step 244). After the 
cascaded processing is complete, the restored signal is delivered to the 
CCM 48 to initiate despreading and demodulation (step 246). Processing of 

20 the next received signal is then initiated. 

FIGS 26-28 are spectrum diagrams illustrating a cascaded 
suppression procedure. FIG. 26 represents a received signal after being 
transfomied to a frequency domain representation. As illustrated, the 
received signal includes both narrow band jammers 250 and frequency 

25 selective fading 252, Adaptive inverse weighting is first applied to remove 
the narrow band jammers 250 and the spectrum of FIG. 27 results. An 
inverse whitening function is then applied to the signal to remove the 
frequency selective fading 252 and the spectrum of FIG. 28 results. Virtually 
any combination of interference suppression methods can be implemented 

30 in the targeted interference suppression unit 78 in accordance with the 
present invention. 
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The principles of tlie present invention can be used in a wide range of 
applications, such as commercial paging, commercial cellular and PCS, 
spread spectrum wireless LAN, satellite uplink/crosslink/downlink, small unit 
operations, digital battlefield, code division multiple access (CDMA), and 
5 others. In one embodiment,, the invention is implemented in full duplex, 
handheld communicators in a mobile communications system. In another 
embodiment, the invention is used in a spread spectrum PCS system which 
is overlaid on existing commercial/military narrowband systems to provide 
simultaneous frequency reuse. In a further embodiment, the invention is 
10 used in a spread spectrum digital paging system which is overlaid on 

existing commercial/military narrowband systems. Other applications are 
also possible. 

It should be appreciated that the interference suppressor 316 of the 
present invention can be used without the interference classifier 314. That 
15 is, criteria other than precise interference type can be used to determine 
which of the interference suppression methods within the interference 
suppressor 316 should be used to perfomi suppression. 
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CLAIMS 

What is claimed is: 

5 1 . An apparatus for reducing interference components in a 

received signal in a spread spectrum communications system, comprising: 

a signal receptor for receiving a signal from a communications 
channel; 

a modal moment estimator for determining modal moment estimates 
1 0 for components in said signal, said modal moment estimates including a 
mean magnitude estimate and a standard deviation estimate; and 

an interference suppressor for suppressing an Interference 
component in said signal using said modal moment estimates. 

15 2. The apparatus, as claimed in claim 1, wherein: 

said modal moment estimator includes means for obtaining a 
frequency domain representation of said signal, said frequency domain 
representation indicating a magnitude level for each of a plurality of 
frequencies in said signal. 

20 

3. The apparatus, as claimed in claim 2, wherein: 

said modal moment estimator includes means for deriving a 
magnitude distribution using said frequency domain representation, said 
magnitude distribution indicating how often each of a plurality of magnitude 
25 levels are reached in said frequency domain representation. 

4. The apparatus, as claimed in claim 3, wherein: 

said modal moment estimator includes means for locating modes 
within said magnitude distribution, each mode corresponding to a signal 
30 component. 
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5. The apparatus, as claimed in claim 1, wherein: 

said interference suppressor includes means for determining an 
adaptive limit value using said modal moment estimates for a first mode. 

6. The apparatus, as claimed In claim 5, wherein: 

said adaptive limit value equals a weighted sum of said mean 
magnitude of said first mode and said standard deviation of said first mode. 

7. The apparatus, as claimed in claim 5, wherein: 

said interference suppressor includes means for applying an inverse 
weighting to first portions of said frequency domain representation that 
exceed said adaptive limit. 

8. The apparatus, as claimed In claim 5, wherein: 

said interference suppressor includes means for excising portions of 
said frequency domain representation that exceed said adaptive limit. 

9. A method for reducing interference components within a 
received signal in a spread spectrum communications system, said method 
comprising the steps of: 

receiving a signal from a channel; 

determining modal moment estimates for components within said 
signal, said modal moment estimates including a mean magnitude estimate 
and a standard deviation estimate; and 

suppressing an interference component in said signal using said 
modal moment estimates. 
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10. The method, as claimed In claim 9, wherein said determining 
step includes: 

obtaining a frequency domain representation of said signal, said 
frequency domain representation indicating a magnitude level for each of a 
5 plurality of frequencies in said signal; 

determining a magnitude distribution for said frequency domain 
representation, said magnitude distribution indicating how often each of a 
plurality of magnitude levels are reached in said frequency domain 
representation; 

0 locating modes within said magnitude distribution, each mode 

corresponding to a signal component; 

finding a mean magnitude for each of said modes; and 
finding a standard deviation for each of said modes. 
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